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Abstract. The normalized incomplete beta function can be defined ei¬ 
ther as cumulative distribution function of beta density or as the Gauss 
hypergeometric function with one of the upper parameters equal to 
unity. Logarithmic concavity/convexity of this function in parameters 
was established by Tinner and Roters in 1997. Their proof is indirect 
and rather difficult; it is based on generalized reproductive property of 
certain more general distributions. These authors remark that these re¬ 
sults ’’seems to be very hard to obtain by usual analytic methods”. In 
the first part of this paper we provide such proof based on standard tools 
of analysis. In the second part we go one step further and investigate 
the sign of generalized Turan determinants formed by shifts of the nor¬ 
malized incomplete beta function. Under some additional restrictions 
we demonstrate that these coefficients are of the same sign. We fur¬ 
ther conjecture that such restrictions can be removed without altering 
the results. Our method of proof also leads to various companion re¬ 
sults which may be of independent interest. In particular, we establish 
linearization formulas and two-sided bounds for the above mentioned 
Turan determinants. Further, we find two combinatorial style identities 
for finite sums which we believe to be new. 
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1. Motivation and introduction. The beta distribution is perhaps the single most 
important continuous compactly supported probability distributions. Its particular cases 
include the Wigner semicircle, the Marchenko-Pastur and the arcsine laws; it is important 
in Bayesian analysis as conjugate prior to binomial and geometric distributions [8]; it plays 

*E-mail: D. Karp - dimkrp@gmail.com 


1 



a role in a large number of applications ranging from population genetics to project man¬ 
agement. The beta distribution is defined by the density — xY~^/B{a,b) supported 

on [0,1], where B{a,b) is Euler’s beta function [3]. The cumulative distribution function 
(CDF) of the beta density is given by the normalized incomplete beta function 

T(„h\ = = s; 

’ B(a,b) 

In a recent paper |1] containing the mathematical analysis supporting a recent investigation 
related to lottery frauds in Florida, the authors solve the following optimization problem 

m m 

CiO^b ^ min under constraint /(ai, as, • • •, = JJ Ip.{wi,ai -Wi + l)>e. 

i i=l 

For efficient numerical solution of such problem it is desirable that the function / be quasi¬ 
concave which implies that the feasible set of the problem is convex and local minimizer is also 
global. Since monotonic transformations do not alter quasi-concavity it suffices then to show 
that log / is a concave function. This fact, in turn, is implied by log-concavity of 6 —)■ /a, (a, b) 
on (0, cxd). The purpose of this paper is to investigate the log-concavity/convexity properties 
of the function /^.(a, 6) viewed as the function of parameters a and b. In particular, we 
demonstrate that b —)■ Ix{a,b) is indeed log-concave on (0,cxo), while a —)■ Ix{a,b) is log- 
convex on (0, oo) for b G (0,1) and log-concave for 6 > 1. Log-concavity of 6 —)■ b) is 

equivalent to the positivity of the generalized Turan determinant 

Ja;(a, b + a)Ix{a, b + P) - ^(a, b)Ix{a, b + a + P) 

for all a, /3 > 0, while log-convexity (log-concavity) of a —>■ /^(a, b) is equivalent to negativity 
(positivity) of 

Ix{a + a, b)Ix{a + I3,b) - Ix{a, b)Ix{a + a + I3,b). 

We go one step further and study the signs of the power series coefficients (in powers of x) of 
the above Turan determinants. Under some additional restrictions we demonstrate that these 
coefficients are of the same sign. We further conjecture that such restrictions can be removed 
without altering the results. Our method of proof also leads to various companion results 
which may be of independent interest. In particular, we establish linearization formulas and 
two-sided bounds for the above Turan determinants. Further, we find two combinatorial 
style identities for hnite sums which we believe to be new. 

When this paper was nearly completed we discovered that log-concavity/convexity of 
the CDF of beta distribution has been demonstrated in 1997 by Finner and Roters in their 
fundamental work . Partial results in this direction have been previously given by the same 
authors in [6] and by Das Gupta and Sarkar in |9]. However, their log-concavity proofs are 
very involved and indirect - they appear as a by-product of generalized reproductive property 
of certain probability measures. Finner and Roters note in their paper that ”we are now 
able to conclude a result for the Beta distribution which seems to be very hard to obtain by 
usual analytic methods”. The proofs presented in this paper are precisely ”by usual analytic 
methods”. The power series coefficients of the Turan determinants have not been considered 
by the above authors, so our results in this direction strengthen the achievements of [7]. 
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We conclude the introduction by presenting an alternative expression for the normalized 
incomplete beta function in terms of the Gauss hypergeometric function [3] 






n =0 

Indeed, using Euler-Pochhammer integral representation 

r(7) 


2 ^ 1 ( 0 ;,/?; 7 ; x) = 


^(1 —^ ^(1 —-ua;) °‘du 


r(/3)r(7 - /3) 7o 

we get by a change of variable and application of Euler’s formula B{a, b) = r(a)r(6)/r(a+6): 

^ - 6, a; a + 1; x). 

f (a + 1)1 ( 0 ) 

Further, applying another Euler’s formula 

2^i(a,/3;7;x) = (1 - 7 -/?; 7; 

we obtain the well-known representation 

{^-xf2Fi{a+b, l;a+l;x) = ^ ^ (1) 


r(a + i)r(6) 


^ r(6)r(a -h 1 n)‘ 

n =0 V / V / 


2 . Log-concavity of Ix{(i, b) in a and b. The proof of the next theorem has been inspired 
by the log-concavity proof for incomplete gamma function given in [1]. See a related result 

in|2|. 

Theorem 1 For each fixed a > 0 and x G (0,1) the function b —)■ Ix{a,b) is strictly 
log-concave on (0, cxd). 

Proof. We need to show that 

r\2 

— log(4(a,6)) < 0. 

It is easy to compute 

^ log(5(a, b)) = fi'{b) - fi'{a -|- 6) > 0, (2) 

where V’(^) = r'( 2 ;)/r(z) is logarithmic derivative of the gamma function and the inequality 
follows from the fact that z is decreasing on (0, 00 ) according to the representation 




~i{n + zfi' 

n =0 ^ ' 


Further, 


02 


= 5a.(a,6) J F ^{l-tf ^[\og{l - t)]‘^dt- 


F ^{l—ty ^\og{l — t)dt 


Uo 


Bx{a, bYififb) - fi'{a b)) =: Ua,b{x) 
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We have t4,fe(0) = Ua^hiX) = 0- The second equality can be seen by setting x = 1 and then 
differentiating which is legitimate by continuous double differentiability or by setting x = 1 
in the above formula and computing the integrals. We aim to demonstrate that Ua,b{x) is 
decreasing on (0,a;i) and increasing on (xi,!) for some Xi G (0,1). This will prove that 
Ua,b{.x) < 0 for a; G (0,1) for each a, 6 > 0. We have 

Va,b{x) := ^ - t)^“^[log(l - t)fdt + B^{a,b)[\og{l - x)]"^ 

X t-*- “ Jo 

— 21 og(l —x) f {1 — \og{l — t)dt — 2Bx{a,b){^jJ'{b) — ^jJ'{a + b)). 

Jo 

It is rather straightforward to see that I4,fe(0) = 0 and 14,6(1) = +oo. We aim to demonstrate 
that Va^b(x) is decreasing on ( 0 , 0 : 2 ) and increasing on (x2, 1) for some x^ G (0,1). This will 
imply that it changes sign on ( 0 , 1 ) exactly once from minus to plus which in turns implies 
the same conclusion for lB^^{x). Taking the next derivative we obtain after cancelations: 

Wa,b(,x) ■= ^(1 - x)Vlf,{x) = -B^{a,b) log(l - x) 

+ [ t^-\l-tf-Hog{l-t)dt-x^-\l-x)\^/J'{b)-^/J\a + b)). 

Jo 

It is clear that II4^f,(l) = +cxd. The value at zero depends on a: if 0 < a < 1 then 
Wa,bi0) = — 00 ; if a = 1 then 114 , 6 ( 0 ) = —(tlj'ib) — 'ijjja + 6 )) < 0 ; if a > 1 then ITa,b( 0 ) = 0 . 
We will demonstrate below that Wa,b{x) is increasing on (0,1) when 0 < a < 1 and Wa,b{x) 
is decreasing on ( 0 , 0 : 3 ) and increasing on ( 0 : 3 ,1) for some 0:3 G (0,1) when a > 1. This will 
imply that it changes sign on ( 0 , 1 ) exactly once from minus to plus which in turns implies 
the same conclusion for Vj^^{x). Taking the next derivative we get: 

Za,b{^) ■= (1 - x)W'^^^{x) = B^{a, b) - ( 4 '( 6 ) - ip'{a + 6 ))o:“"^(l - o:)^((a - 1)(1 - 0 :) - bx). 

We have Za^[,(l) = B(a, b) > 0. At x = 0 five cases reveal themselves: 

Case I: if 0 < a < 1 then Za^/,(0) = +cxd. In this case Za^i,{x) > 0 on (0,1) since 
((a — 1)(1 — 0 :) — bx) < 0 and ip'{b) — ip' {a + 6) > 0; 

Case II: if a = 1 then Za,b{0) = b{ip'{b) — ip'{a + b)) > 0; We will show below that 
Z'afkx) = 0 on (0,1) so that Za,bix) = b{ip'{b) — ip'{a + b)) > 0 for all x G (0,1); 

Case III: if 1 < a < 2 then Za,b{t)) = — 00 ; 

Case IV: if a = 2 then Za,b{0) = —ijJ'ib) — ip'{a + b)) < 0; 

Case V: if a > 2 then Za,f,(0) = 0; 

Taking one more derivative we get: 

Qa,b{^) ■■= ^',fe(o:)o:^"“(l - 0:)^"'’ 

= x^ — {ip'{b) — Ip'{a + b)){{a + b — l)‘^x‘^ — (a — l)(2a + 2b — 3)x + (a — l)(a — 2)). 
Straightforward calculation yields: 

Qa,6(0) = -(pj'ib) - ip'{a + b)){a - l)(a - 2), Qa,b{^) = 1 “ b'^{p)'{b) - ip'{a + b)). 


4 



Since a —)■ Qa,b(l) is decreasing on (0, cxd), Qo,b(l) = 1 and Qi,b(l) = 0 because '^'(l + b) = 
'ip'{b) — l/&^, these formulas lead to the following conclusions. 

Case I; if 0 < a < 1 then Qa,b(0) < 0 and Qa,fe(l) > 0 for all 6 > 0. Since Qa,b{x) is a 
quadratic this implies that it has exactly one change of sign on ( 0 , 1 ). 

Case II: if a = 1 then Qa,bix) = 0 on [0,1] for all 6 > 0 since ijj^l + b) = ijj'ib) — 1/6^. 

Case III: if 1 < a < 2 then Qa,fe(0) > 0 and Qa,b{^) < 0 for all 6 > 0. Since Qa,b{x) is a 
quadratic this implies that it has exactly one change of sign on ( 0 , 1 ). 

Case IV: If a = 2 then Qa,fe(0) = 0 and Qa,b{^) < 0 for all 6 > 0. It follows that 
Qa,b{x) > 0 on (0, a) and Qa,b{,x) < 0 on (a, 1) for some a G (0,1). An apparent alternative 
Qa,b{x) < 0 on ( 0 , 1 ) cannot hold since Zafi{x) would then be decreasing on ( 0 , 1 ) while 
Za,bi0) < 0 and ^a, 6 (l) > 0 rules out such possibility; 

Case V: a > 2 then Qa,b(0) < 0 and Qa,b{^) < 0 for all 6 > 0. This implies that Qa,b{x) 
follows the sign pattern on (0,1)- An apparent alternative Qa,b(x) < 0 on (0,1) 

cannot hold since Za^/,(x) would then be decreasing on ( 0 , 1 ) while Za^h(O) = 0 and 2 ’^ ^( 1 ) > 0 
rules out such possibility; 

Looking carefully at each case we see that for all 6 > 0 the function Za,b{x) is positive on 
( 0 , 1 ) ifO < a < 1 and changes sign from minus to plus if a > 1. This implies that lL((^(x) 
behaves in the same way, so that in all cases Wa,b{x) changes sign exactly once from minus 
to plus. This in turn implies that V^f^{x) has the same behavior and hence Vafi{x) is hrst 
decreasing and then increasing. In view of the boundary values at x = 0 and x = 1 this 
means that 14,b(x) changes sign exactly once from minus to plus. Finally, it follows that 
U'^jj{x) is negative on (0,xi) and positive on (xi, 1) for some Xi G (0,1), yielding Uafi{x) < 0 
on ( 0 , 1 ) as claimed. □ 

Dehnition of /^(a, b) immediately leads to the reflection formula 

4(a,6) = 1 - Ji_„(6,a). 


This implies 

— log(4(a,6)) = [4(a,6)]“^pi_,,(6,a)—/i_a;(6,a) - 

\Ix{,Oj-iby\ ^^2 n). 

The hrst term is negative according to Theorem [TJ The second term, however, may change 
sign depending on the values of b and x as demonstrated by numerical evidence. Hence, we 
cannot draw any dehnitive conclusion about log-concavity of a —)■ /a,(a, b) from the rehection 
formula. Indeed, it turns out that the result depends on the value of b. Nevertheless, a 
method similar to that used in the proof of Theorem [1] also works here combined with the 
following lemma. 

Lemma 1 For fixed a > 0 define /q(x) : (0, oo) — )■ (— oo, cxd) by 

fa{x) = {i^ix + a) - 'iIj{x)Y + fi'{x + a) - fi'{x). 

Then fa{x) <0 if 0<a<l and /a(x) >0 ifa>l. 


-^Ii-xib,a) 

da 
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Proof. Using recurrence relations '0(x + 1) = 'ip{x) + 1/x and 'ip{x + 1) = ijj{x) — 1/x^, we 
calculate 


f{x + l) - f{x) = 

(V’(a; + Q! +1) —'0(a; +1))^ — {'i/j{x + a) — -ilj{x)y + -ip' {x + a + 1) —ij^x + a) +i/j'{x) —i/j^x + I) 
= {tlj{x + a + 1) — + 1) + + a) — V’(a^))('0(a^ + a + 1) — tfj^x + 1) — + a) + ^J^x)) 


x’^ {x + a )2 


= ( 2i){x + a) + 


X + a 


2'ip{x) - 




XJ \x + a XJ x‘^ (x + a)2 


a 


x{x + a) 


2'ip{x) — 2'ip{x + a) + 


a 


a{2x + a) 
x{x + a) J x‘^{x + ay 
2a 


'ipix) — ipix + a) H— . 

X 


x{x + a) 

Writing ga{x) for the function in parentheses and utilizing the integral representations 

1 


tjj^x) = -7 + 




1 — e 


— t 


dt and 


X 




where 7 stands for Euler-Mascheroni constant, we get: 


daix) = 


.00 g-xt^g-at _ g-t 


1 — e 


—t 


-dt. 


The last formula makes it obvious that 5 fQ,(x) > 0 for 0 < a < 1 and ga{x) < 0 for a > 1. 
Further, from the asymptotic formulas 

ij^x) =log{x) -+ 0{x~^), tl^^x) = - + ^ + 0{x~^), X ^ 00, 

2x X 2x^ 

we conclude that lim f{x) = 0. Altogether this implies for 0 < a < 1: 

X^OO 

f{x) < f{x + 1 ) < f{x + 2 ) < • • • < lim f{x + m) = 0, 

m—)-oo 

and for a > 1 : 

f{x) > f{x + 1 ) > f{x + 2 ) > • • • > lim f{x + m) = 0 , 


which proves the lemma. □ 

Stronger results for a function similar to fa{x) but still slightly different from it can be 
found in [T71 Theorem 1 . 2 ], 

Theorem 2 Suppose x G (0,1) is fixed. Then: 

(a) f /0 < 6 < 1 the function a —)■ Ix{a, h) is strictly log-convex on ( 0 , 00 ). 

(b) ifb = l, the function Ix{a,b) = and so is log-neutral. 

(c) if b > 1 the function a —?■ /^(a, b) is strictly log-concave on ( 0 , cx)). 
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Proof. The case 6=1 can be verified directly. It remains to show that 

— log(/a;(a, 6)) > 0 if 0 < 6 < 1 and — log(4(a, 6)) < 0 if 6 > 1. 

For completeness, however, we will inclnde 6=1 into the forgoing considerations. By 
symmetry we have from ([2]) 

^ log(5(a, 6)) = - ^l)'{a + 6) > 0. 

Farther, 

B^{a,hY—\og{I^{a,h)) = B^{a,h) j - tf~^[\og{t)fdt- 

^“-1(1 _ log(t)dt 

Direct verification yields Ua^iS^) = Ua^iX) = 0. Next, we have 


- B^a, bf{X{a) - + 6 )) =: Ua,b{X- 



KbX) 


VaXX ■= a-in -^ ^(1 - tf ^[\og{tXdt + 5,,(a,6)[log(a:)] 

X Jo 


— 21og(a;) / 6“ ^(1 — t)” ^ \og{t)dt — 2BXa,b){X{X ~ XX + b)). 


Repeatedly using L’Hopital’s rule we compute the limits: 

BXa,b) 1 (log(l/a:))^ 

hm -r — , ^ = hm ——;—-—r = - hm- 

a;-:>o [log(a:)] 2 2(log(l/a:)) ^ 2 x “ 


= 1 3(log(lA))^ ^ 3 21og(lA) _ 3 1^0. 


2 a :^0 ax “ 


2 x^o a^x “ 


x^o a'^x 


3 /y> — CL 


Similarly, 


lim log(a;) / 6 “ ^(1 — 6 )* ^\og{t)dt = 0, 
x^o L 


so that I4,fe(0) = 0. Further, 


Rl - ^[log(6)]^(i6 = i?(a, 6)[('0(a) -'0(a + 6))^ +'^'(a) -+ 6)], 


so that 

I4,fe(l) = B{a,b)[XX) -' 4 ’{a + b)f - (XX) -XX + X)]- 

It follows from Lemma [D that (a) 14,6(1) < 0 for 0 < 6 < 1; (b) I4,b(l) = 0 for 6 = 1; (c) 
14,6(1) > 0 for 6 > 1; We aim to demonstrate that 14,^(a;) has precisely one change of sign 
in cases (a) and (c) and is identically zero in case (b). Taking the next derivative we obtain 
after cancelations: 


X, 


Wa,bX) ■= -jV^bX) = 5„(a,6)log(a;) 


^a-l(i _ ^^6-1 _ xf-^{XX) “ XX + X). 
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The boundary values are: Wa,b{0) = 0 and 


Wa,b{l) 


— CXD, 0 < & < 1 

0 , b = l, 

B{a, b){ip{a + b) — ip{a)) >0, b > 1. 


These values follow from the evaluation 

^“- 1(1 _ iog(t)(it = B{a, b){ilj{a) — ijj{a + b)) 

and the recurrence relations + 1) — t/’(a) = 1/a, + 1) — = —1/a^. Taking the 

next derivative we get 



Za,bix) = xW^ ,,{x) := Bx{a, b) + (i/j'ia) - + 6))a;“(l - x)'" ^((6 - l)x - a(l - x)). 

Clearly, Za^h(0) =0. At a; = 1 hve cases reveal themselves: 

Case I: if 0 < 6 < 1 then ^a,b(l) = —oo; 

Case II: if 6 = 1 then ^a,fe(l) = 0; We will show below that Z'^i^(x) = 0 on (0,1) so that 
Za,b{x) = 0 for all x E (0,1); 

Case III: if 1 < 6 < 2 then ^a,b(l) = +oo; 

Case IV: if 6 = 2 then ^a,6(l) = B{a, b) + ijj'{a) — 'ip^a + 6) > 0; 

Case V: if 6 > 2 then ^a,b(l) = B{a, b) > 0. 

Finally, computing one more derivative and making some rearrangements we arrive at 


Qa,b{x') 




X' 


a—1 


(1 -x) 


b-3 


= (1 — xf — (ip'ia) — 'ip\a + b)){{a + b — — (( 2 a + 1)(6 — 1 ) + 2a?)x + or) 


The boundary values are 

Qa,fe( 0 ) = 1 - o?{'il)'{a) - Ip'{a + 6 )), Qa,b{^) = -(^'(a) - ^'(a + b)){b - l){b - 2 ). 

Since b —)■ (5a,b(0) is decreasing on (0, cxd), (5a,o(0) = 1 and (5a,i(0) = 0 because ■^'(1 + a) = 
'ip'i^a) — 1 /a^, these formulas lead to the following conclusions. 

Case I: if 0 < 6 < 1 then Qa,b{0) > 0 and Qa,b{^) < 0 for all a > 0. Since Qa,b{x) is a 
quadratic this implies that it has exactly one change of sign on (0,1). 

Case II: if 6 = 1 then Qa,b{x) = 0 on [0,1] for all a > 0 since ■^'(1 + a) = 'ip'i^a) — 1/a^. 

Case III: if 1 < 6 < 2 then (5a,fe(0) < 0 and Qa,bi^) > 0 for all 6 > 0. Since Qa,bix) is a 
quadratic this implies that it has exactly one change of sign on (0,1). 

Case IV: If 6 = 2 then (5a,fe(0) < 0 and Qa,bi^) = 0 for all a > 0. It follows that 
Qa,b{x) < 0 on (0, a) and Qa,b{x) > 0 on (a, 1) for some a E (0, 1). An apparent alternative 
Qa,b{x) < 0 on (0,1) cannot hold since Za,b{x) would then be decreasing on (0,1) while 
Za,b{0) = 0 and Za,b{^) > 0 rules out such possibility. 

Case V: b > 2 then (5a,6(0) < 0 and Qa,b{^) < 0 for all a > 0. This implies that Qa,b{x) 
follows the sign pattern (—,+,—) on (0,1). An apparent alternative Qa,b(x) < 0 on (0,1) 
cannot hold since Za,b{,x) would then be decreasing on ( 0 , 1 ) while Za,b{0) = 0 and Za,b{^) > 0 
rules out such possibility. 



Tracing back the primitives Qa,b —^ ^a,6 —^ ^ K,6 ^ Ua,b we see that for all a > 0 

the function Za,b{x) changes sign from plus to minus on (0,1) if 0 < 6 < 1 is identically zero 
if 6 = 1 and changes sign from minus to plus if 6 > 1. This implies that Wa,b{x) and 14,b(x) 
behave in the same way, so that Ua,h{x) is hrst increasing and then decreasing for 0 < 6 < 1, 
is identically zero if 6 = 1 and is hrst decreasing and then increasing if 6 > 1. In view of the 
boundary values at a: = 0 and a: = 1 this means that Uafi{x) > 0 on (0,1) for 0 < 6 < 1 and 
Ua,b{x) < 0 for 6 > 1 as claimed. □ 

3. Coefficient-wise log-concavity of Ix{a,b). As mentioned in the introduction, loga¬ 
rithmic concavity of 6 —)■ Ix{a,b) demonstrated in Theorem [T] is equivalent to the inequality 

(f>a,b,a,^{x) := Ix{a, b a)Ix{a, b + /3) - Ix{a, b)Ix{a, b + a + /3) > 0 (3) 

valid for all positive values of a, b, a and (3 and each x G (0,1). This inequality expresses 
the strict Wright log-concavity of 6 —)■ Ix{a,b) - a property equivalent to log-concavity for 
all continuous functions [HI Chapter 1.4], [T^ Section 1.1]. On the other hand, formula (fTj) 
implies that the function of (j)a,b,a,/ 3 {x) defined in ([3]) has the power series expansion 

OO 

0a,6,«,/3(a:) = a;^“(l - ^ (t)kX^. (4) 

k=0 

It is then natural to ask whether and when the coefficients (pk are positive. In general, for a 
given formal power series 

OO 

f{f^',x) = '^fk{n)x’' (5) 

fc =0 

with nonnegative coefficients /^(/i) which depend continuously on a parameter fi from a 
real interval /, we will say that fi —)■ f{fi]x) is coefficient-wise Wright log-concave on I for 
a G A C [0, cxo), G i? C [0, cxd) if the formal power series for the product difference 

f{fi + a; x)f{fi + l3-,x)- /(/r; x)f{fi + a + l3;x) 

has nonnegative coefficients at all powers of x for a & A, f3 ^ B and all fi E I. If coefficients 
are strictly positive we say that the corresponding property holds strictly. If coefficients 
are non-positive we talk about coefficient-wise Wright log-convexity. These concepts using 
slightly different terminology were introduced by the author (jointly with S.I. Kalmykov) in 

ng. 

Conclusions of Theorem [2] can be restated in terms of the function 

OO 

'ipa,b,a,p{x) := Ix{a + a,b)Ix{a + /3,b)-Ix{a,b)Ix{a + a + l3,b) = 3 ;^^+"+/^{1 -xf^ y^^kx’" (6) 

k=0 

as follows: if a > 0 and 0 < a; < 1 then 'ilJa,b,a,i 3 {x) < 0 for 0 < 6 < 1 and 'ilia,b,a,i 3 {x) > 0 for 
6 > 1. Just like with (j)a,b,a,i 3 {x) it is then natural to consider the coefficient-wise log-convexity 
and log-concavity of a —)■ — x)~^Ix{a, b), i.e. study the sign of the coefficients ipk- 

The main purpose of this section is to give partial solution to the following conjectures 
representing substantial strengthenings of Theorems [T] and [2], respectively. 
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Conjecture 1 The coefficients (l)k are positive for all a,h,a,(3 > 0 and x G (0,1), so 
that h —)■ — x)~^Ix{a,b) is strictly coefficient-wise Wright log-concave on (0, cx)) for 

these values of parameters. 

Conjecture 2 The coefficients ipk are negative for all a, a,/? > 0 and x G (0,1) z/0 < 
h <1 and positive ifb > 1, so that a —>■ x~°‘{l — x)~^I^{a,b) is strictly coefficient-wise Wright 
log-convex on (0, oo) for a, /3 > 0, x G (0,1) and 0 < b < 1; and a —)■ — x)~^Ix{a, b) is 

strictly coefficient-wise Wright log-concave on (0, oo) for a,/9 > 0, x G (0,1) and b > 1. 

Our approach also leads to a number of related results which may be of independent 
interest. These include alternative formulas for the coefficients fk and ipk, linearization 
identities for the product differences of hypergeometric functions, two-sided inequalities for 
the Turan determinants formed by normalized incomplete beta functions and two presumably 
new combinatorial identities. 

We will need the following lemma (see m Lemma 3]). 

Lemma 2 Suppose p. —)• /(p;x) is coefficient-wise Wright log-concave on [0, oo) for 
a = 1, (3 > 0, i.e. 

f{p -1- 1; x)/(/i + /3; x) - /(p; x)/(/i + /3 + 1; x) 

has nonnegative coefficients at all powers of x for all pffi > 0. Then p —)■ f{p;x) is 
coefficient-wise Wright log-concave on [0, oo) for aGN, /5>a — 1. 

Theorem 3 Suppose a G M. Then the following identity holds true 

a,b,affix) _ r(a + h + j)r(a + b-\-l3-\-a — j — 1 )^^ 

^ 2 a+a+y^l_^yb [T{b)]W{a + j + l)T{a + l3 + a-j) 

( a b j _ \ a-\-b-\-(3-\-a — 1—j l'l,a-\-b-\-(3-\-a—j _ 

\ a -|- j -|- + j -\-2 ’y a -\- /3 -\- a — j ^ + l3-\-a—j-\-l ’ yj 


Proof. First we investigate the function 'ipa,b,iffiffi dehned in ([H]). Using ([T]) we compute 


ffi,bpffiffi [r(&)]^ _ T{a + b)T{a + b + p) f a + b f l,a-^b-^l 


x2“+/^+i(l — x)2^ r(a -I- l)r(a + p-\- 1 ) a -I- + 2 


2F1 


l,a-\-b-\- p 
a p 1 


a b -\- p a -\- b 

-2x1 


, PI l,a-{-b-\-p-{-l 


) ^ j r ^ ^ X 

' ' m=0 


X 


E 

k=0 


{a -\- b)(a -\- b -\- l)fc(a p)m-k {a -\- b -\- p)(a -\- b)k{a S b -\- p -\- l)m-fc 
(a -|- l)(a -|- 2)k{a S p -\- l)m-k (a -|- /i -|- l)(a -|- l)fc(ci + p -\- 2)m,-k 

_ — 1) (® + b)k{a + b -\- p)m-k 




{a + l){a + p + 1) (a + 2)k{a + p + 2)m-k 


{m — 2k + p) 


= A, 




d -f- b /ijd-hb-l-l 


Cl “h 1 


>Fi 


a -|- 2 


;x - 


a -\- b -\- p 

a -\- p -\-1 \a p 2 


;x 


X 
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where we have used the summation formula 


m 


E 


(a + b)k{a + b + fi)m-k 
{a + 2)fc(a + /i + 2)m-k 


{m—2k+fi) 


(ci + l)(a + /r + 1) 

i^) 


(a + b)m+i 
{a + l)m+l 


{a + b + fi)m+i 
(a + /i + l)m+l 


demonstrated in [m p.72]. Next we set g{a) = Ix{a,b)/x°'{l — x)^ and calculate 


'4^a,b,a,/3 


— x) 


2b 


g{a + a)g{a + (3) — g{a)g{a + a + (3) 


=[g{a+a)g{a+l3)-g{a + a-l)g{a+l3 + l)] + [g{a+a-l)g{a+l3 + l)-g{a+a-2)g{a+l3 + 2)] 
+• ■ ■-\-\g{pi-\-2')g{(i-\- [3-\-oi —2) — g{(i-\-V)g{(i-\-j3-\-oi —l)] + [5'(u+l)(y'(u+/3+Q;—1) — g{a3)g{a-\-j3-\-(y.y\ 
_ ^ A I g + 6 + j /l,a + 6 + j + l 

[r(fe)]" ^ a+i,M,«+/3-2i-l I a + + J+2 

a + b-\-(3-\-Oi — 1—j /l,a + 6 + /3 + a — j _ \1 

a-\-(3-\-OL — j ^ ^yg + Z^ + a —j + 1 ’ j j 



Substituting the value of Aa+j^b,i,a+/ 3 - 2 j-i into this formula yields ([7]). 

The next corollary refines Theorem [2] and gives partial solution to Conjecture |2J 


□ 


Corollary 1 Suppose a G N. Then the coefficients fjk defined in ([6]) ean be computed by 


! \ r(g + 6 + j)r(g + 6 + /3 + o — j — 1 ) 

^ ^ [r(&)]2r(g + 3 + i)r(g + /? + «- j) 

,, f (g + + j)k+i (g + & + /3 + Q — 1— j)k+i ) /gN 

\(g + j + l)fc+i {a + (3 + a — j)k+i J 

Furthermore, if 0 < b < 1 then fjk < 0 for all a,/3 > 0 and if b > 1 then ipk > t) for 

/d>a — 1 , /3>0 and a G N. 


Proof. Formula dH]) follows immediately from ([7]). The case 0 < 6 < 1 falls under conditions 
of [ini Theorem 3] and the claim follows from this theorem. Further, according to Lemma [2] 
for /3 > a — 1 the sign of fjk = 'ipk{c(,(3) is determined by the sign of 'ipk{3,f3)- In view of 
(g)fc = F(g + fc)/F(g) for a = 1 we have 

, F(g + 6)F(g -\- b -\~ jo) F(g -\- k -\- 6)F(g + 1) F(g -\- ja -\~ k -\- 6)F(g + /i + 1) 

^ ^ ’ F(g + l)F(g + yU + 1) F(g + 6)F(g + A; + 1) T{a + p, + b)T{a + pt, + k + 1) 

Define the function 

, T{a + p + k + x)T{a + x) 

F(g + k + a;)F(g + p + x) 

This function is decreasing according to [5l Theorem 6]. It is straightforward that fjk-ii^, (3) < 
0 is equivalent to /(I) < f{b) while fjk-ii^, (3) > 0 is equivalent to /(I) > fib), so that the 
claim follows from the decrease of f{x). □ 

This corollary can also be derived from our previous results in uni Theorem 3] and 
m Theorem 2], however, the above derivation has the advantage of presenting the explicit 
formulas for the coefficients tjjk which do not follow from these references. 
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Remark. Substituting ([T]) for into the definition of 'il’a,b,a ,/3 given in (|6]) and changing 
notation we can rewrite ([7]) as the following linearization identity for the quadratic form in 
Gauss hypergeometric functions 


T\ -2-^1 

(c)n+l 


o + n + l,! \ _/a + /3,l 

c + n + 1 + 


j 


X 


_ {a + (3)r,+i fa, 1 _ fa + /? + n + 1,1 

(c + /3Ui^ V ^^)^^\c + p + n + l 




X 


_ (Q)j(Q + /3)n-j f 0- + j p/a + j + 1,1 \ 

+ \c + f \c + J + l ’ J 


a + l3 + n- j ^ a + P + n + 1 - j,l 
c + P + n-j'^ \c + (3 + n + l—j 



(9) 


Here n = a — 1 is any nonnegative integer, while other parameters as well as x can be 
arbitrary complex by analytic continuation. In particular when n = 0 the above identity 
reduces to (after some rearrangement) 


a 

—2 

c 


Fi 


1, Cl + 1 
c + 1 



2 FI 


1, a + /? 

C + /3 



a +13 fl,a + (3 + l 

7Tp^^\c + f^ + i 





1 


Formula ([9]) leads to a curious combinatorial identity which is presumably new. 


Corollary 2 For arbitrary nonnegative integers m, n the following identity holds 


(a)fc(Q + (d)m-k r (g + k)n+i _ {a + (d + m — k)n+i | 

“ if>)kip + ld)m-k I (& + k)n+i {h + (d + m- k)n+i J 

_ (g)j(g + (d)n-j r (g + j)m+l _ {a + fd + n — j)m+l 1 

\ib + j)m+i (b + (d + n — j)m+i ) 

Proof. The claimed identity follows on equating coefficients at x^ on both sides of the (p|) 
and using the formula (p)„+i(/i + n + 1)^ = {n)k{F + ^)n+i- □ 

The next corollary presents two-sided bounds for the Turan determinant formed by shifts 
of /a;(g, b) in the first parameter. 

Corollary 3 Suppose 6 > 1, z/ G N and a > 12 — 1. Then 

mx°‘{l — xY < Ix{a, hY — Ix{a — u, b)Ix{a + I'yb) < MIx{a, bY (10) 

for all X G (0,1). Here 

[r(a -|- b)Y r(a + b — i')T{a + b + u) 

"" ^ [r(6)]2[r(a + l)]2 “ [r(6)]2r(a-z/ + l)r(a + z/ + l)’ 

^ _ (g l)Ya + b - - {a + 1 - v)Ya + b)^ 

(a + l)^(a + b- 
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Proof. Inequality from below is just a rewriting of the inequality 'il^a-v,b,u,v > 

which follows from nonnegativity of the coefficients 'ipk for all A; = 0,1,... To prove the upper 

bound consider the function 


a —)■ G{a) 


I^{a,b)T{a + 1) 
x“(l — x)^r(a + b) 


1 (q + b)n n 

r(6 )^(a + l). 


According to [121 Theorem 2] the above function is log-convex for 6 > 1, so that for a > v — 1 
we have G{a — h')G{a + i') > [G(a)]^. This inequality can be rewritten as 


Ix{a,b- u)Ix{a,b + u) > 


[r(a + l)]^r(a + b + u)r{a + b-iy) ^ 
[r(a + 6)]2r(a + z/ + l)r(a -u + l) ^ 


Simple rearrangement of the above inequality gives the upper bound in (ITU]) . □ 

To investigate the power series coefficients of the function (f>a,b,a,i 3 {x) defined in (|3D we 
need the following lemma. 


Lemma 3 The following identity holds: 


+ b)k{a + b + I3)m-k ^ qI , , + ^ +/^)m-|-l ~ (« + ^)m+l) 

- . , [o[ZK — m) + p{a -t- K}} — - —- 

^_Q T T l)m—A: T l)m 


Proof. Denote 


Uk = ^ +^)m k ^ fc = 0,1,..., m. 

(a -|- l)fc(a -|- l)m-k 

An application of Gosper’s algorithm yields anti-differences: 

a^(a -|- 6)fc(a (3)m-k+i 


'^k — 


, k = 0,1,... ,m + 1. 


(®)fc (®)m—fc+l 

Straightforward substitution provides verification of this formula. Hence, we get 

E a((a -I- 6 -|- (3)m+i — (a + fe)m+i) , 

Uk 'Uq ^ • n 


/c=0 


(a -|- l)r 


Theorem 4 Suppose a G M. Then the following identity holds true 


4^a,b,a,P (2^) 


ol—1 


ar(a + b + j)r(a + b-\-f 3 -\-a—j —1) 
2,2a(1 _ ^yb+a+y [r(a + l)]2r(6 + j + l)r(6 + ^ + a-j) 


X 


(^a -\- b (3 a — j — 1)2^1 


l,a + b + (3 + a- j 

fl -|- 1 


;x - (a + b + j) 2 Fi 


l,a-|-6 -|-jTl 
Cl -|- 1 


■,x 


( 11 ) 
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Proof. First we investigate the function defined in ([3]). We have 


0a,6,1,(3 ( 2 ^) 


2.2a (1 _ 2;)2fe+l+/3 


00 m ^ 

m=0 k=0 ^ 


r(a + 6 + 1 + k)T(a + b + (3 + m — k) 
i r(6 + l)r(a + 1 + k)T{b + 0)r(a + 1 + m - k) 


r(a + 6 + fc)r(a + b + (3 + l + m — k) 
r(6)r(a + 1 + k)T{b + 1 + /9)r(a + 1 + m — k) 


= 

m=0 


X 


r(a + b + k)T{a + b + /3 + m — k) 

^ r( 6 )r (6 + 0 )r(a + l + fc)r(a + l + m-k) 


E 


fc =0 


a + b + k a + b + f3 + m — k 


T{a + b)T{a + b + P) 


(g + &)fc(a + & + (3)m-k 

r (6 + l)r (6 + 0 + l)r(a + 1)2 ^ ia + l)kia + l)m-k 


E^^E 


6 + /? 


(6(2A; —m)+/3(a + /c)). 


Application of Lemma [3] then gives 

0a,fe,i,/3 gr(a + 6 )r(a + 6 + /3) 


X 


2a 


(1 - a;) 2 ^+i +/3 r(b + i)r (6 + 0 + i)r(a + 1)2 ^ 


(g + 6 + /3)m+i — (g + &)m+i 


m=0 


gr(g + 6 )r(g + 6 + 0 ) 
r (6 + i)r (6 + 0 + i)r(g + i )2 


)g + 6 + 0)2^1 


(g A l)m 

l,g + 6 + 0 + l 
g + 1 


X ‘ 


X 


-(a + 6 ) 2^1 


1 , g + 6 + 1 
g + 1 
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Next denote /( 6 ) = /a;(g, 6 )/x“(l — x)^. For a G N we have 


+ “)/('> + / 3 ) - /( 6)/(6 + o + / 3 ) = l/(f> + a )/(6 + /?) 

— /(6 + a — l)/(6 + 0 + l)] + [/(6 + a — l)/(6 + 0 + l)—/(6 + a — 2)/(6 + 0 + 2)] + -- - 
+ [/(^+2)/(6+0+a—2)—/(6+1)/(6+0 +q;—1)] + [/(6+1)/(6+0+q;— 1)—/(6)/(6+0+a)]. 

Each expression in brackets has the form x~'^'^(l—x)~'^^~‘^~^(/)a^b^.a-j,i,/ 3 -a+ 2 j-i, j = 1, 2,..., a. 
Hence, we can apply formula to each such bracket which yields flTT]) . □ 

The next corollary refines Theorem [1] and gives partial solution to Conjecture [1] 

Corollary 4 Suppose a G N. Then the coefficients fk defined in (jl]) can be computed by 


gF(g + 6 + j)F(g + 6 + 0 + a — j — 1) 

^ ^ [r(a + 1)]2F(6 + J + 1)F(6 + 0 + a - J)(g + 1)/ 

{(g+ 6 + 0 + 0 — j — l)fc+i — (g + 6 + j)k+i} ■ (13) 


If, furthermore, 0 > a — 1 and 0 > 0 then 0^ > 0. 


Proof. Formula fll3p follows immediately from fllip . According to Lemma H] for 0 > o — 1 
the sign of 0^ = 0fc(a, 0) is determined by the sign of 0 a:( 1, 0). In view of (IT^ for o = 1 we 
have 


0 fc(l, 0 ) 


gF(g + 6 )F(g + 6 + 0) 

F (6 + 1)F(6 + 0 + l)F(g + l)^(g + l)fc 


((g + 6 + 0 )fc+i 


(g + 6 )fc+i), 
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which is clearly positive for /3 > 0 . □ 

Remark. Substituting ([T]) for into the definition of (()a,b,a,i 3 given in (|3]) and changing 
notation we can rewrite (|TT |1 as the following linearization identity for the quadratic form in 
Gauss hypergeometric functions 


(ci + &)n+i „/ci + fe + n + l,l \ /a + 6 + /3,1 \ 

—nr\ -2^^! II ; ^ 2i^i ,1 ; X 

{b)n+i \a + l J \a + l ) 

(ci + 6 + (5')n+i + 1 \ _/o + 6 + /3 + 'n + l,l 

■ (t+/ju. ^ v+i -y^^v+i 

a, {a, -\- h)ji^a, -\- h -\- /3)n—j 

+ /^) j^Q + l)j(^ + /^ + l)n-i 


■|^(ci + 6 + 


P + n- j)2Fi 


a + b + /3 + n + l— j,l 

Cl “h 1 



(a + b + j) 2 Fi 


a + 6 + j + l,l 
Cl “h 1 



(14) 


Here n = a — 1 is any nonnegative integer, while other parameters as well as x can be 
arbitrary complex by analytic continuation. In particular, when n = 0 the above identity 
after some rearrangement reduces to 


-rT 2 ^i 

a — c + 1 \c 



■x] 1 1 

[ V 

’7 6-c+lJ 


b - c+1 


iFi 


1,6 + 1 


:x 


2 F 1 


1 , a 


;x + 


c — 1 

a — c + 1 


Formula ffT411 leads to a curious combinatorial identity which is presumably new. 
Corollary 5 For arbitrary nonnegative integers m, n the following identity holds 


a 

bib + /3) 


E 


i=o 


(a + 6)j(a + 6 + f3)n-j 

{b + l)j{b + 1 + (3)n-j 


{{a + b + (d + n 


j)m+i — (a + 6 + j)m+l} 


— (a + l)ra 

k=0 


(a + 6)fc(n + 6 + (3)m-k 

{a + l)fc(a + l)m-k 


(a + 6 + k)n-\-i 
(6)n+l 


{a + b + /3 + m- k)n+i 


{b + (d) 


n+l 


Proof. The claimed identity follows on equating coefficients at x^ on both sides of the flTT)) 
and using the formula (/r + n + l)fc(/i)n+i = (+)fc(+ + k)n+i- □ 

The next corollary presents two-sided bounds for the Turan determinant formed by shifts 
of Ix{ci, b) in the second parameter. 


Corollary 6 Suppose u is a positive integer b > u, a > 0. Then 

ma:^“(l - < [4(a, h)f - 4(a, b + b - v) < M[4(a, b)f (15) 

for all X G (0,1). Here 

[r(a + b)Y r(a + 6 — v)V{a + b + v) 

"" ^ [r(a + i)]2[r(6)]2 “ [r(a + i)]2r(6-i/)r(6 + z/)’ 

_ {b)v{a + b-v)y-{a + b)^{b - v)^ 

{b)y{a + b-v)^ 


15 






















Proof. Inequality from below is just a rewriting of the inequality (j)a,b-Ly,u,u > 

which follows from nonnegativity of the coefficients (pk for all A; = 0,1,... To prove the upper 

bound consider the function 


b F{b) 


Ua,b)T{b) 
x“(l — x)^r(a + b) 


(a + b)n (a + b)nn\ x^ 

-“Tfa + n + l)^ “r(a + n + l)n! 

n=0 ^ ^ 71=0 ^ ' 


It is easy to check that the coefficients {n!/r(a + n + l)}j^o form a log-convex sequence. 
Then according to [121 Theorem 1] the above function is log-convex, so that for b > u we 
have F{b — v)F{b + v) > [F{b)Y. This inequality can be rewritten as 


Ix{,a,b- v)lYa,b + p) > 


[T{b)YT{a + b + iy)T{a + b-iy) 
[r{a + b)Yr{b + u)r{b - u) 


[Ix{a,b)Y. 


Simple rearrangement of the above inequality gives the upper bound in ([TB]) . □ 
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